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Abstract. We report the detection of high excitation CO emission from the most distant quasar currently known, 
SDSS J114816.64+525150.3 (hereafter J1148+5251), at a redshift z = 6.419. The CO (J =6->5) and (J =7~*6) 
lines were detected using the IRAM Plateau de Bure interferometer, showing a width of ~ 280kms~^. An upper 
flux limit for the CO (J =1— >0) line was obtained from observations with the Effelsberg 100-meter telescope. 
Assuming no gravitational magnification, we estimate a molecular gas mass of « 2 X 10^° Mq. Using the CO 
(3^2) observations by Walter et al. (2003), a comparison of the line flux ratios with predictions from a large 
velocity gradient model suggests that the gas is likely of high excitation, at densities ~ 10^ cm~^ and a temperature 
^ 100 K. Since in this case the CO lines appear to have moderate optical depths, the gas must be extended over 
a few kpc. The gas mass detected in J1148-I-5251 can fuel star formation at the rate implied by the far-infrared 
luminosity for less than 10 million years, a time comparable to the dynamical time of the region. The gas must 
therefore be replenished quickly, and metal and dust enrichment must occur fast. The strong dust emission and 
massive, dense gas reservoir at z ~ 6.4 provide further evidence that vigorous star formation is co-eval with the 
rapid growth of massive black holes at these early epochs of the Universe. 

Key words, galaxies: formation - galaxies: starburst - galaxies: high-redshift - quasars: emission lines - quasars: 
individual: SDSS J1148-I-5251 - cosmology: observations 



1. Introduction et al. 2003). Optical, radio and millimeter observations 

?H indicate that Jl 148-1-5251 could be weakly amplified by 

; The luminous quasars at redshifts z > 6 found in the Sloan intervening lens (Fan et al. 2003; White et al. 2003; 

Digital Sky Survey by Fan et al. (2001, 2003) provide a g^^^^j^. 2003), but in what follows, we will assume 

unique opportunity to study the formation of massive ob- ^^^^ ampHfication. J1148+5251 is an very luminous 

jects during the epoch at which the intergalactic medium ^^^^^ ^^^^^^ ^ 2^ 3^ ^^^^ ^ ^^m^^^ p^^^^^^ ^ 

was being reionized by the first luminous sources (Becker g^pgrniassive (« 3 x 10^ Mq) black hole radiating close 

et al. 2001; Kogut et al. 2003; Cen 2003). Studying signa- -^^ Eddington luminosity (WiUott et al. 2003). If the 

tures of star formation in these exceptional objects is also ^^^^ ^^^^ ^^^^^^ ^^^^ associated with J1148+5251 

of great mterest to test whether the correlation between . „™„^^.;„„„i 4.„ >,i„„u >,„ip „„„„ ;„ „ similar 

Id [jl (j[JKJL LiUlidi LU hlLfZj UldCiV llUlc illcLoD 111 d W dV blllllldl 

the central black hole mass and the stellar bulge mass ob- ^^^^ ^^^^^ .^^ ^^^^^ spheroids (Shields et al. 2003), 

served in local spheroids (Magorrian et al. 1998; Gebhardt -^^ ^^^^ ^^^^^ he ^ 2 x lO^^ Mq, and J1148+5251 would 

et al. 2000) can be traced to the early formation stages ^^^^^ ^^^^ ^^^^-^^ collapsed structures to have 

of quasars and their host galaxies. ^^^^^^ ^^^1 Universe (e.g., Haiman & Loeb 2001). 

J1148-h5251, at a redshift oi z = 6.42 (Fan et al. 
2003), is the most distant quasar known, observed only 

« 850 million years after the Big Bang (we adopt Ho = detection of thermal dust emission in 

71 km s-i Mpc-i, flA = 0.73 and fl^ = 0.27 - Spergel J1148+5251 (Bertoldi et al. 2003) imphcs a far-mfrared lu- 

minosity of « 10 i© and a dust mass of 7 x 10 Mq. If 

* Based on observations obtained with the IRAM Plateau de the dominant heating mechanism is radiation from young 

Bure Interferometer, and with the Effelsberg 100 m telescope, stars, then the star formation rate implied from the FIR 
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luminosity is ~ 3000 A/q yr^^ which requires vast amounts 
of molecular gas to be maintained. 

Molecular gas in excess of 10^" Mq was detected 
through their CO emission in fifteen z > 2 far-infrared ul- 
traluminous (Lfir > 10^^ Lq) radio galaxies and quasars 
(e.g., Cox, Omont, & Bertoldi 2002). At z > 4, CO emis- 
sion was detected towards four quasars (Omont et al. 1996; 
Ohta et al. 1996; Guilloteau et al. 1997, 1999; Cox et 
al. 2002). The CO emission was resolved in BR 1202-0725 
(CariUi et al. 2002) at z = 4.69, the highest redshift CO 
detected so far, and PSS 2322-^1944 at z = 4.12 (Carilh 
et al. 2003). The extended nature of the CO provides the 
most direct evidence for active star formation in the host 
galaxies of distant quasars, and indicates that black hole 
accretion and star-formation are closely related. 

To explore the growth of massive black holes and 
their associated stellar populations at the end of the 
"dark ages", we have searched for CO emission toward 
J1148-I-5251. We here report the detection of CO (6^5) 
and (7^6) line emission. In a separate study, Walter et 
al. (2003) report the discovery of CO (3— ^>2) emission us- 
ing the Very Large Array (VLA). 

2. Observations 

Observations of the CO (7^6) and (6— >5) emission lines 
were made with the IRAM Plateau de Bure interferome- 
ter between March and May 2003. We used the 6 antenna 
D configuration which results in a beam of 5.7" x 4.1" at 
3.2 mm. The 3 mm receivers were tuned in single sideband 
and the typical SSB system temperatures were w 150 K. 
The total integration time was 14 hours for CO (6— >5), 
and 22 hours for CO (7— >6). The spectra are displayed 
in Fig. ^ and the image of the averaged data is shown 
in Fig. 13 Within the astrometric uncertainties of ±0"3, 
the CO emission coincides with the optical position given 
by Fan et al. (2003). At the 5" resolution of our Plateau 
de Bure observations, the CO emission is unresolved, con- 
sistent with the VLA C0(3— >2) detection, which is unre- 
solved at 1'.'5 (Waher et al 2003). 

The CO (6^5) and (7^6) lines are detected at 
centroid frequencies corresponding to a redshift 6.419 
(Table . Within the uncertainties this agrees with the 
z = 6.41 ±0.01 of the Mgii A2799 line (Willott et al. 2003) 
(Fig.^. The CO redshift, which is likely to correspond to 
the systemic redshift of the quasar, differs significantly 
from the range z — 6.36 — 6.39 derived from high ion- 
ization UV lines (White et al. 2003), which trace high 
velocity (> 1000 kms^^), blue-shifted gas related to the 
quasar activity. 

No continuum emission was detected in our coadded 
3 mm data, which includes observations at other frequen- 
cies, but excludes the continuum redward of the CO lines, 
where weak line emission may be present. At the position 
of J1148-I-5251 we obtain a continuum flux of 0.09 ± 0.13 
mjy. At 43 GHz the continuum remains undetected with 
—31 ± 57 /iJy. These upper limits are consistent with the 
measured 250 GHz flux density of 5.0 ± 0.6 mJy (Bertoldi 
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Fig. 1. J1148-I-5251 spectra of CO (6^5), (7^6), and 
their average, binned to 64, 55, and 55 kms~^, respec- 
tively, four times the original spectral resolution. Zero ve- 
locity corresponds to the centroid of the (6^5) line at 
93.206 GHz. Gaussian fits with FWHM= 279km s^^ are 
shown as light lines. 
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Fig. 2. Velocity-integrated (±220kms~'^) map of the av- 
eraged CO (6-^5) and (7^6) emission. Contour steps are 
0.34 mJy/beam = 2(7. The cross indicates the optical po- 
sition. 



et al. 2003), if we adopt a grey body spectrum with tem- 
perature > 50 K and dust emissivity cx v^. 

We fit Gaussians to the line spectra within ±300 
kms~^ of the centroid, with no baseline subtraction. The 
best fit line widths of the three spectra shown in Fig. 1 
range between 280 and 320 kms^^, similar to the widths 
found in other high redshift quasars (e.g., Cox et al. 2002). 
The width and centroid of the best Gaussian fit to the 
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Table 1. Properties of the CO lines observed toward SDSS J1148+5251 
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[mJy] 


[kms-1] 


[Jy km s""*"] 


[10"'Kkms-ipc^] 


[10*^0] 


CO (7-^6) 


806.651 


108.729 


6.4192±0.0009 


2.14 


279I" 


0.64±0.088 


1.73±0.24 


2.92±0.40 


CO (6^5) 


691.473 


93.204 


6.4187±0.0006 


2.45 


279 


0.73±0.076 


2.69±0.24 


2.86±0.25 


CO (3-*2) 


345.795 


46.610 


6.419±0.004 


0.6 


320* 


0.18±0.02 


2.68±0.27 


0.35±0.04 


CO (1-+0) 


115.271 


15.537 




< 0.36 




< 0.11+ 


< 14.2 


< 0.070 



NOTE. - For J1148+5251, the apparent CO line luminosity is given by L'qq — 3.2 x 10'*/coJ^obs' intrinsic line luminosity 
Leo = 4.2 X 10^/co !/obs, in the units given above (see Solomon et al. 1997). Upper hmits are 3a. + Adopting the line width of 
CO (6-^5). * Line width corresponds to the 50 MHz channel width of the VLA 46.6 GHz observations. 



(6^5) line were adopted for the Gaussian fit to the lower 
quality (7^6) line to determine its flux. The CO (6^5) 
and (7^6) line fluxes are determined at strong confldence 
levels of 10(7 and 7a, respectively. 

We searched for CO (1^0) emission using the 
Effelsberg 100-meter telescope in March and April 2003 
with a 1.9 cm HEMT receiver (T^ys - 40 K on a 
scale, aperture efficiency ~40%, beam width 60", position 
switching mode). The integration time was ~50 hours, 
yielding a r.m.s. of 0.4 mK (~ 0.4 mJy per 24 km 

s~^). No CO (1^0) emission was detected at the redshift 
found for the higher CO transitions (Table 

3. Discussion 

To constrain the physical conditions of the molecular gas 
in J1148+5251 we compared the observed CO line flux 
ratios with those predicted by a one-component large ve- 
locity gradient (LVG) model (Mao et al. 2000). The fine 
flux ratios are determined by the gas density, temperature, 
and the optical depth in the CO lines, i.e., the column 
density of CO per velocity interval. The large flux ratio 
between the (6^5) and (3^2) lines implies that the gas 
has a high excitation. The lower excitation of the J = 7 
level suggests a moderate optical depth (Fig. 3). 

High gas densities are typical of the molecular gas 
present in the nuclear regions of nearby starburst galax- 
ies (Solomon & Downes 1998). The most extreme con- 
ditions so far were found in the starburst nucleus of 
NGC 253 (Fig. 3), where the CO excitation is similar to 
that in J1148+5251. Detailed LVG modeling by Bradford 
et al. (2003) indicate that the CO, ^^CO, and H2 data of 
NGC 253 are consistent with n(H2) = 4.5 x lO^cm"^ ^nd 
T = 120 K, and CO line optical depths r < 4 (Fig. ^. 
With only three line fluxes and one upper limit, we cannot 
constrain the physical conditions of the gas in J1148-f5251 
as tightly: the data can be fit both with low-opacity mod- 
els, in which temperature and density are degenerate, 
Xn}/'^ Ki 2.5 x lO"' K cm^^-^, and with high-opacity lower- 
excitation models. With the high excitation temperatures 
for the gas in NGC 253 and J1148-I-5251 the cosmic back- 
ground temperature (3 K and 20 K, respectively) does not 
affect the gas excitation notably. 

To infer the total molecular mass from the CO emis- 
sion in the absence of constraints on the CO abundance, 
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Fig. 3. Integrated line flux, /cO: normalized to CO 
(6— >5). Diamonds show the values for J1148-I-5251. The 
dashed line shows line flux increasing as , which is 
expected for optically thick conditions. The solid lines 
show LVG models with Tki„ 120 K and ^(Ha) = 4.5 x 
lO^cm"'^, with different maximum optical depth in the 
CO lines. The dotted line shows the line flux distribution 
observed for the starburst nucleus of NGC 253 (Bradford 
et al. 2003). 

one typically adopts an empirical conversion factor, a, be- 
tween the apparent CO (1^0) line luminosity L'qq, and 
the total molecular mass, Mjja ■ However, this conversion 
may depend on the CO excitation. In nearby starbursts 
with moderate gas excitation, Downes & Solomon (1998) 
derive a = 0.8 M0 (K kms~^ pc^)~^. The high excita- 
tion and moderate line opacities of our one-component 
LVG model for J1148-f5251 predict a CO (l->0) line 
flux much lower than that of an optically thick distri- 
bution (Fig. 3). As observed for NGC 253, it is hkely 
that an optically thick, low-excitation molecular compo- 
nent adds to the lower J level populations. Given this un- 
certainty, and since a value of a for very high excitation 
conditions is unknown, we do not estimate the mass on 



4 



F. Bertoldi et al.: CO in a quasar host at z = 6.42 



the one component LVG predictions, but extrapolate to 
-^co(i^o) ^ '^■'^ ^ 10^°Kkms~^pc^, with the assumption 
of a constant Hne brightness temperature (the optically 
thick case) from J = 1 to 6. With the quoted conversion 
factor, we find A/hj ~ 2 x 10^" M©. We then estimate the 
gas to dust mass ratio Mnj/Mdust ~ 30, which is similar 
to the values found for local ULIRGs and other high red- 
shift quasars (e.g., Guilloteau et al. 1999; Cox et al. 2002). 

The minimum area of the molecular region can be 
estimated from the ratio of the observed line brightness 
temperature (11 niK for CO (6— >5)) and the intrinsic 
line brightness, which is 23 and 56 K in the LVG mod- 
els with T = 4 and 1, respectively (Fig. With a 5" 
beam the corresponding source radius (assuming uniform 
coverage) is 0.1-0.15", or 560-840 pc. Placing 2 x IO^^Mq 
in a volume of radius 560 pc gives an average CO col- 
umn density 1.2 x 10^°(Xco/10~'*)cm~^, where Xco is 
the CO abundance relative to H2. If the observed CO 
lines have moderate optical depth, r < 4, the CO col- 
umn is < 2 X lO^^cm"^, which would either require a 
low CO abundance, Xco < 2 x 10~^, or a larger volume. 
Considering the low gas-to-dust ratio estimated above and 
the high metalicities implied by the optical lines (Fan et 
al. 2003), a CO abundance much lower than the Galactic 
~ 10~^ seems unlikely. Rather, a larger radius of ^ 1400 
pc for the gas distribution could account for the moderate 
CO line opacities. 

If the molecular gas forms an inclined disk (angle i rela- 
tive to the sky plane) in Keplerian rotation about a spher- 
ical mass, the line width and a minimum source radius 
between 560 and 1400 pc yield a minimum gravitational 
mass enclosed by the disk of (2 — 6) x 10^ sin~^ i Mq. For 
large inclination angles this mass would not be much larger 
than that of the black hole, and a factor 4-10 smaller than 
the gas mass implied by the line intensities. The latter may 
have been overestimated given the approximate nature of 
our estimate; alternatively, the CO disk inclination is close 
to the sky plane, i ~ 20 — 30 deg, which is more likely con- 
sidering the large dust mass, and the fact that the AGN 
is optically unobscured. 

The detection of large amounts of dense molecular gas 
in J1148+5251 supports the conjecture that the strong far- 
infrared luminosity seen from many quasars arises from 
extended star forming regions, and is not due to heat- 
ing from the AGN (Omont ct al. 2001, 2003; Carilh et 
al. 2001). Although for Jl 148-1-5251 the emission remains 
spatially unresolved, the large masses of warm CO are un- 
likely to be heated by the AGN at a kpc distance. With 
the estimated mass of molecular gas of ^ 2 x 10^" Mq 
star formation in J1148-I-5251 could be sustained at the 
rate ~ 3OOOM0yr^^ implied by the far-infrared luminos- 
ity for a short time only, < 10 million years. This is 
comparable to the estimated duty cycle time of quasars 
(e.g. Wyithe & Loeb 2003), and to the dynamical time of 
the star forming region, which implies a rapid gas deple- 
tion unless the system continues to accrete gas at a high 
rate. If the replenishing gas is of low-metalicity, the short 



depletion time suggests that the enrichment with heavy 
elements and dust is rapid, which leaves only supernovae 
and winds from the most massive stars as possible sources. 

Our low estimate for the dynamical to luminous mass 
ratio excludes the presence of a large stellar mass within 
the volume of the CO emission. The duration of star for- 
mation at the present rate could therefore not have been 
much longer than 10^ yr, unless the starburst does not 
form many long-lived, low-mass stars, in which case the 
star formation rate would have been overestimated and 
the depletion time could be longer. 

The dynamical mass is an order of magnitude smaller 
than the bulge mass deduced from the correlation between 
the black hole mass and the bulge mass (or velocity dis- 
persion) in local spheroids (Magorrian et al. 1998). This 
could be due to a biased selection of a non-representative, 
bright quasar, or it confirms a tendency for the stellar to 
black hole mass ratio to decrease at higher redshifts (Rix 
et al. 1999), possibly due to self-regulating star formation 
mechanisms (Wyithe & Loeb 2003). 
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